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Abstract 


This  report  discusses  the  Cloud  Rise  Module  (CRM)  of  the  Defense  Land  Fallout  Interpretive 
Code  (DBLFIC),  the  DOD’s  reference  fallout  code.  The  first  section  discusses  the  errors  found  in 
the  1979  documentation.  These  errors  have  been  corrected  and  further  research  by  the  Air  Force 
Institute  of  Technolog)'  will  include  the  corrections  listed  in  this  section.  The  second  section  will 
discuss  the  use  of  Mathematics  in  simulating  the  CRM.  This  part  of  the  report  will  describe  how  a 
liigher  level  language  was  used  to  review  tlie  results  of  DELFlC’s  CRM  and  uncover  some  prob¬ 
lem  areas.  The  final  part  of  this  report  will  describe  further  work  as  planned  by  the  author  in 
improving  the  CRM  of  DELFIC. 
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Critique  of  DELFIC’s  Cloud  Rise  Module 

1.0  Introduction 

This  report  is  t  critique  of  the  Defense  Land  Fallout  Inteipretive  Code  (DBLFIC)  as  it  per¬ 
tains  to  the  author’s  research  into  nuclear  cloud  rise  and  growth.  The  version  of  DELFIC 
being  reviewed  in  this  report  is  the  last  well  documented  change  which  gave  reasonable 
re.sults  (Norment,  1979a  tt  1979b).  Other  efforts  have  gone  into  DELFIC  since  then  but 
have  either  not  been  well  documented  or  give  questionable  results.  The  1979  version  was 
recovered  by  McGahan  with  SAIC  and  provided  to  AFTT  fcr  this  research.  Throughout 
this  report,  “the  documentation”  refers  to  Volume  I  (Norment,  1979a)  and  Volume  II  (Nor¬ 
ment,  1979b)  of  the  1979  version  of  DELFIC. 

The  current  AFTT  research  is  focused  on  the  Cloud  Rise  Module  (CRM)  of  DELFIC.  The 
CRM  includes  the  main  subprogram  ICRMEX  and  its  subordinate  routines,  which  deter¬ 
mine  the  initial  conditions,  solve  the  cloud  rise  equations  and  prepare  the  definition  of  the 
particles  aloft  for  the  diffusive  transport  module  (DTM).  The  code  in  the  CHIM  was  com¬ 
pared  to  the  listings  in  Volume  n. 

This  report  has  three  parts.  The  first  section  discusses  the  errors  found  in  the  documenta¬ 
tion.  These  errors  have  been  collected  and  further  research  by  AFTT  will  include  the  cor¬ 
rections  listed  in  this  section.  The  second  section  will  discuss  the  use  of  Mathematics  in 
simulating  the  CRM.  This  part  of  the  report  will  desoibe  how  a  higher  level  language  was 
used  to  review  the  results  of  DELFIC’s  CRM  and  uncover  some  problem  areas.  The  final 
part  of  this  repon  will  describe  further  work  as  planned  by  the  author  in  improving  the 
CRM  of  DELFIC. 

2.0  Corrections  Needed  in  the  1979  DELFIC  CRM 

The  last  well  documented  version  of  DELFIC  is  Norment’s  1979  version.  In  studying  this 
version,  a  few  errors  were  found  both  in  the  volume  pertaining  to  the  theory  and  the  vol¬ 
ume  which  included  the  source  code.  Below  is  a  discussion  of  those  discrepancies  and  the 
corrections  needed  to  remedy  them.  When  reviewing  these  errors,  it  may  be  beneficial  to 
have  a  copy  of  Norment’s  two  volumes  from  1979  on  hand. 
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2.1  Changes  needed  to  compile  CRM 

In  recovering  the  1979  version  of  DELHC,  McGihan  had  to  add  a  MAIN  program,  since 
the  documentation  was  lacking  one.  Since  AFTT’s  research  is  focused  on  the  CRM  of 
DELFIC,  calls  in  the  MAIN  program  to  the  DTM  or  (he  output  processor  module  (0PM) 
were  commented  out.  The  CRM  compiled  with  no  changes  on  a  386-SX>16  using 
Microsoft  FORTRAN  5.1.  To  compile  on  a  SPARC  station  2  using  SUN  FORTRAN  1.4, 
two  separate  FORMAT  statements  (016  and  099)  in  the  DBG  subroutine,  had  to  be  edited. 
The  line  numbers  as  listed  in  Volume  n  were  DBG  25-DBG  27  and  DBG  30-DBG  34.  To 
prevent  ambiguity  in  the  statements,  commas  were  needed  after  the  X’s,  which  represent 
spaces  in  FORTRAN  output.  These  changes  allowed  the  CRM  to  be  compiled  without  any 
warnings  or  errors. 

Similar  changes  were  needed  to  compile  the  entire  DELPIC  program  but  will  not  be  listed 
here.  The  documentation  contains  a  test  case's  input  and  output  which  was  used  to  verify 
that  the  code  was  running  correctly.  After  the  CRM  was  compiled  successfully,  the  pro¬ 
gram's  ou^ut  was  verified  with  the  output  as  listed  in  Volum'*  n.  Additional  post-proces¬ 
sors  were  built  with  Mathematics  to  display  the  CRM  variables  as  a  function  of  time.  The 
code  is  fast  on  today's  computers  with  the  entire  DELFIC  test  case  taking  only  20  seconds 
on  a  486-33. 

2.2  Errors  in  1979  Volume  I 

There  were  two  errors  in  Volume  I.  Both  pertained  to  the  definition  of  the  symbol  for  the 
ratio  of  molecular  weights  of  water  to  air,  The  errors  are  described  in  the  following  sec¬ 
tions  with  the  heading  giving  the  page  number  in  Volume  I. 

2.2.1  Definit'on  of  ^  pages  14  and  93 

^  is  incorrectly  defined  as  the  ratio  of  molecular  weights  of  air  to  water,  29/18.  In  fact, 
after  checking  with  earlier  documentation  as  well  as  other  references,  the  constant  ^  is  that 
of  water  over  air.  18/29.  This  is  only  an  error  in  Volume  I  and  therefore  the  code  has  the 
correct  definition  for 

2.2.2  Definition  of  x,  page  13 

The  above  error  affects  the  definition  of  x^,  Eq.  (2.1.11). 
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(EQI) 


_  HrP., 

t  §/> 

The  ^  in  this  equation  should  be  in  the  numerator  not  the  denominator.  This  is  the  only 
occurrence  of  misusing  ^  since  this  is  the  only  new  equation  in  the  1979  version  that  uses 
All  the  other  equations  which  use  ^  were  derived  in  earlier  versions  (Heubsch,  1967; 
Norment,  1970;  Norment,  1977)  with  its  proper  definition.  These  other  equations  are 
repeated  in  their  correct  form  in  the  1979  version. 

This  error  was  found  using  Mathematica.  3y  coding  the  equations  as  they  were  presented 
in  Volume  I,  it  was  noticed  that  the  value  for  Xe  was  a  factor  ^  too  low.  This  caused  more 
detailed  scrutiny  of  the  equation  until  the  error  was  found. 

23  Errors  in  1979  Volume  II 

The  errors  in  Volume  n  have  to  do  with  the  code  listings  for  the  CRM.  The  errors  are 
described  in  the  following  sections  with  the  heading  giving  the  line  number  of  the  subrou- 
tine  which  contains  the  error  along  with  the  page  number  in  Volume  n.  All  corrections 
were  made  individually  to  check  the  effect  on  the  1979  test  case  (a  50  kt,  0  ft.  HOB  shot). 
The  results  are  shown  along  with  the  original,  uncorrected  results  (see  Figure  ORIGI¬ 
NAL).  h  should  be  noted  that  even  when  the  effect  on  this  test  case  is  negligible,  the  effect 
of  the  error  may  be  greater  with  other  scenarios. 

23.1  CRMIN66page78 

7  soilht=3sam*(tad+78 1 .6*(tpr-te)+0.2856*(tpr**2-te**2)+ 

The  documentation  gives  the  defimtions  for  the  specific  heat  of  soil,  Cj,  on  page  13  of  Vol¬ 
ume  I.  When  the  code  does  the  integration  of  the  definition  for  Cj,  the  term  0.5612  T 
should  go  to  0.2806  T^.  Instead  the  code  lists  it  as  0.2856  as  shown  above.  All  that  needs 
to  be  done  is  to  change  the  5  into  a  0.  This  has  little  effect  on  the  test  case. 

23.2  DERIV  91  page  85 

qq=qt*qx*qxe*(  1  .+x+wt)/(  1  .+x+s+wt) 

In  the  code  P’  is  incorrectly  calculated  (l.+X+WT)/(l.+X+S+WT)  as  shown  above.  The 
numerator  should  be  (l.+X),  consistent  with  the  true  definition  for  p’(sce  page  92  Volume 
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I).  The  CRM  uses  this  vtriable  in  defining  the  first  term  on  the  right  hand  side  of  Eq. 
(2.2.4)  on  page  21  of  Volume  I. 


1  ^  _g 

2mdt„~  mdt 


This  change  had  little  effect  on  the  test  case. 


(EQ2) 


233  DERIVU3page86 

100  dnn=(nn/(l.-cpai/(cp*t*qx)))*rmix*(rs  •rl+(qt*qx*qxe*9.8*u-eps)* 

This  line  of  code  is  used  to  calculate  the  rate  of  mass  entrained  as  defined  in  Bq.  (2.2  JD) 
on  page  22  of  Volume  I. 


dm 

tat 


P' 


P'm 


-(pHv+ 


P7 

‘  t 


P'  ft*  1 


(EQ3) 


The  error  occurs  in  the  fraction  before  the  integral  sign.  The  code  as  listed  above  deletes 
the  P’  and  uses  Cp  instead  of  Cp.  To  remedy  this,  the  line  was  replaced  with  the  following 
line,  noting  that  it  is  now  longer  than  the  72  characters  allowed  in  FORTRAN,  requiring 
the  excess  characters  to  be  placed  on  line  DERIV114. 


100  drm=(rm/(l.-rmix*cp8i/(cr*t*qx)))*rmix*(rs  *rl+(qt*q!t*qxe*9.8*u-eps)* 

This  change  has  a  small  but  noticeable  effect  on  the  test  case  (see  figure  DERTV  113  Com¬ 
parison). 


2J.4  ATMR  187  thru  ATMR  203  page  72 
do  260i=2,256 
alt(i)=alt(i-l>Klalt 
225  if(al  .ge.alt(i))go  to  250 
if(alt(i)-al  .It.  2.)  go  to  250 
na=na+l 

if(nat-na  .ge.0)go  to  240 
230  irror=-230 
go  to  130 

240  read(irise)al,a2,a3,a4,a5,a6 
go  to  225 

250  terp=dalt/(al-alt(i-l)) 

8tp(i)=atp(i-l)+terp*(a2-atp(i-l)) 

prs(i)=prs(i-l)+tetp*(a3-prs(i-l)) 

rlh(i)=rlh(i-l)+terp*(a4-ilh(i-l)) 
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FIGURE  2.  Effect  of  error  in  (otraiudMit  equatioa,  DERIV  U3  CompariMW 


rho(i)=rho(i-l)+teip*(»5-rhc(i-I)) 

eto(i)=eu(i-l)+terp*(»6-et»(i-l)) 

260  continue 

The  above  lot^  expands  the  atmosphere  Ubles  to  256  entries  from  -1,000m  to  SO, 000m  in 
200m  increments.  It  is  supposed  to  be  using  linear  interpolation  of  the  input  data  to  do 
tlus,  but  it  isn't  In  reality  it  is  linear  inteipolating  between  the  last  200m  increment  value 
it  has  calculated  and  stored  and  the  next  input  datum.  This  causes  a  skewing  of  the  data  in 
the  table,  which  is  used  quite  extensively  throughout  the  QIM. 

The  best  way  to  correct  this  is  to  not  make  an  artificial  regular  spaced  table  to  interpolate 
from,  but  rather,  linear  interpolate  on  the  input  data  itself  when  needed.  This  has  a  small 
effect  on  the  test  case. 


2J.5  DERIV  54  thru  DERIV  61  page  85 
rmix=(  1  .+x)/(  1  .+x+s+wt) 
cr=q)*rmix 
if(tmps-t)380.38 1,381 

381  if(t-848.)3810,3810.3811 

3810  cs=781.6+0.5612»t-1.881e7/t*»2 
go  to  3812 

3811  cs=1003.8+0.13510*t 

3812  cr=cr4cs*(s+wt)/(l.+x+s+wt) 

The  definition  for  Cp  at  the  top  of  page  21  Volume  I  doesn’t  quite  match  the  lines  from  Vol¬ 
ume  n  above.  The  above  code  actually  has  a  switch  in  the  definition  of  ^  which  is  not 
mentioned  in  Volume  I. 


Cp  =  P'Cp+  (1-P')c,t{7',r,)  (EQ4) 

If  the  temperature  of  the  cloud  is  greater  than  the  initial  temperature  of  the  soil  in  the 
cloud,  k(T,  TjIsO.  When  the  cloud  temperature  drops  below  the  initial  soil  temperature, 
k(T,  Tj)=l.  This  is  equivalent  to  stating  that  the  heat  content  of  the  soil  is  negligible  com¬ 
pared  to  that  of  the  air  mixture  until  the  cloud  has  cooled  to  the  initial  soil  temperature. 

If  one  models  the  equations  as  described  in  Volume  I,  always  accounting  for  the  soil  heat, 
a  noticeable  difference  is  seen  in  the  behavior  of  the  cloud  (see  Figure  UERIV  54  Compar¬ 
ison).  This  is  exactly  how  this  difrerence  was  found,  by  modeling  the  Volume  I  equations, 
as  is,  in  Mathematica.  The  difference  seen  using  just  the  test  case  from  1979  shows  that 
ignoring  the  heat  content  of  the  soil  is  not  a  very  good  approximation. 


Crilkpie  of  D£UrCt  aoud  RIm  MoSuk 


II 


OftKiiK  ^TDELFIOi  Ooiid  HUM  MoAili 


3.0  Lessons  Learned  Using  Mathematica 

Mathematica  u  a  higher  level  programming  language  with  symbolic,  numeric,  and 
graphic  capabilities.  By  combining  all  three  formats  in  one  package,  along  with  hundreds 
of  built  in  functions/rouUnes,  this  type  of  language  can  be  used  as  a  separate  debugging 
tool  when  verifying  the  wjrkings  of  a  large  nested,  conventional  program.  It  was  espe¬ 
cially  helpful  since  DELFIC  was  written  by  many  different  authors  using  different  ver¬ 
sions  of  FORTRAN.  Rather  than  trying  to  follow  someone  else’s  programming  from  start 
to  finish,  a  higher  level  language  can  be  used  to  model  the  theory  and  check  for  accuracy. 

When  taking  such  an  approach,  errors  may  be  found  a  little  easier  since  whole  sections  of 
conventional  codes  are  replaced  by  a  few  lines  of  higher  level  code  or  one  built  in  func¬ 
tion  This  approach  also  provides  a  way  of  seeing  better  ways  of  solving  a  problem.  This 
section  describes  some  of  these  findings  not  previously  mentioned  in  the  above  chapter. 

3.1  Method  of  solving  the  set  of  cloud  equations 

One  area  of  DELFIC  in  which  the  higher  level  language  provided  some  insight  was  the 
method  of  solving  the  system  of  coupled  cloud  equations  (see  Appendix  A).  The  set  of 
eight  first  order  nonlinear  differential  equations  are  coupled,  so  that  some  equations  use 
the  derivative  of  another  variable.  DELFIC  uses  a  fourth  order  Runge-Kutta-Gill  method 
with  a  fixed  time  step  for  each  of  three  time  domains.  For  the  first  second,  the  equations 
use  al/32  second  time  step,  then  a  1/4  second  tin.r  step  until  1(X)  seconds,  and  finally  a  2.S 
second  time  step  until  stabilization.  When  transitioning  from  no  co.ndensed  water  to  con¬ 
densed  water  present,  a  1/4  second  time  step  is  used  to  better  define  when  this  transition 
takes  place. 

The  built  in  function  fcr  Mathematica  which  solves  systems  of  DDEs  uses  a  variable  time 
step  and  adapts  to  a  different  method  depending  on  the  stiffness  of  the  equations.  Because 
the  equations  in  DELFIC  are  stiff,  the  Runge-Kutta-Gill  method  may  not  be  a  good  choice 
for  solving  them.  This  may  be  one  area  where  DELFIC  needs  improvement.  The  question 
also  comes  up  as  to  whether  or  not  DELFIC  uses  the  proper  time  steps  for  numerical  accu¬ 
racy,  and  if  so,  is  it  the  best  choice  for  computational  efficiency. 

3.2  IVansition  from  dry  to  wet  equations 

As  mentioned  above,  DELFIC  finds  where  the  onset  of  condensation  occurs  in  the  cloud 
using  a  2.5  second  time  step,  but  then  backs  up  one  step  and  finds  the  onset  more  pre'';:>eiy. 
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It  does  this  by  going  back  to  a  1/4  second  time  step  until  the  transition  and  returning  to  a 
2.5  second  time  step  once  its  found.  This  was  simulated  with  the  nigher  level  language, 
but  a  discrepancy  still  existed  between  the  ou^ut  of  the  two  different  languages. 


It  was  determined  that  the  difference  was  due  to  the  fact  that  DELFIC  fixes  the  equations 
as  either  wet  or  dry  for  the  entire  time  step.  Since  its  a  four  step  method  for  each  time  step, 
the  equations  should  be  allowed  to  switch  as  necessary  within  the  time  step.  This  is  pre¬ 
cisely  what  Mathematica  does.  DELFIC  may  be  trying  to  compensate  for  the  fixing  of  the 
dry/wet  mode  dunng  a  four  step  method  with  its  retrace  using  a  smaller  time  step,  'o  more 
accurately  define  when  the  transition  takes  place. 

3J  Fallout  and  its  effect  on  cloud  rise 

The  way  the  CRM  determines  the  mass  change  due  to  fallout  during  cloud  rise  is  another 
area  of  concern.  In  the  test  case  as  cuirently  modeled  in  the  CRM,  fallout’s  contribution  to 
the  buoyancy  in  the  cloud  is  negligible.  However,  the  documentation  states  that  it  can  be 
important,  and  therefore  the  CRM  spends  a  fair  portion  of  the  computational  effort  to  cal¬ 
culate  it  during  the  solving  of  the  ODEs.  How  it  models  this  was  looked  at  quite  closely 
when  trying  to  emulate  the  model  with  Mathemadca. 

DELFIC  treats  the  cloud  as  if  it  were  a  cylinder  when  it  comes  to  determining  fallout  dur¬ 
ing  the  rise.  For  each  time  step,  the  CRM  determines  what  distance  the  panicles  in  the 
cloud  would  drop.  It  then  reduces  the  soil  mass  by  the  ffacdon  of  the  cylinder  (equiva¬ 
lently  the  fraction  of  the  vertical  height)  which  fell  out  the  bottom  of  me  cloud.  One  must 
keep  in  mir.J  that  DELFIC  assumes  the  entire  mass  is  present  at  its  initialization  tune, 
which  is  a  gross  approximation  to  the  actual  time  dependent  loading  of  the  cloud  with  soil. 
Therefore  allowing  for  a  time  dependent  soil  loss  during  the  rise  seems  very  aniticial  con¬ 
sidering  the  above  assumption.  This  area  will  be  looked  at  more  closely  in  the  future. 

3.4  Oscillating  cloud  height 

DELFIC’s  CRM  doesn't  allow  the  cloud  center  to  drop  from  the  maximum  height 
obtained  and  fixes  its  height  upon  reaching  this  maximum  (see  lines  DERIV  70  to  DERIV 
80  in  Volume  II).  It  was  seen  using  the  Mathematica  model  that  the  cloud  oscillates  in  alti¬ 
tude  before  stabilizing  at  an  altitude  lower  than  the  maximum,  which  DELFIC  freezes  at. 
McGahan  stated  that  oscillation  in  cloud  height  is  a  physical  phenomena  observed  during 
the  atmospheric  tests.  This  was  one  of  the  changes  made  to  the  1979  DELFIC  by  McGa¬ 
han  which  will  be  retained  during  the  research  at  AFIT. 
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3.5  Wind  shear  and  its  elTect  on  cloud  rise 


In  coding  the  equations  in  the  higher  level  language  the  first  time,  the  effect  of  wind  shear 
was  inadvertently  left  out.  This  was  due  to  'he  fact  that  it  was  mentioned  as  a  replacement 
rule  on  page  26,  Eq.  (2.2. 13)  of  Volume  I. 

pliv=»p(£v+^v,)  (EQ5) 

By  comparing  the  results  before  and  after  wind  shear  terms  were  inserted,  it  was  seen  that 
wind  shear  has  a  notic<"<bIe  effect  on  the  test  case  (see  Figure  WIND  SHEAR  Compari¬ 
son).  DELFIC  simply  cai'jolates  the  wind  shear,  Vj,  between  the  veiocity  vector  at  the  top 
of  the  cloud  and  the  bottom.  The  question  remains  if  the  CRM’s  method  of  calculating  the 
wind  shear  is  a  good  approximation  considering  the  vertical  extent  of  some  clouds. 

3.6  Fits  to  physical  data 

A  question  or  two  came  up  when  modeling  the  viscosity  and  specific  heat  equations  used 
in  the  CRM.  Are  they  valid  for  the  temperatures  experienced  in  the  cloud  at  ea'ly  times 
and  are  the  jumps  in  the  specific  heat  fits  valid?  The  viscosity  equation  was  found  to  match 
existing  data  up  to  2000  K  but  data  for  higher  temperatures  was  not  available.  Although 
the  source  of  the  fits  is  still  not  known,  it  was  discovered  that  the  model  is  very  sensitive  to 
the  specific  heat  values.  Equations  for  p,  and  P  will  also  be  checked  because  of  their 
importance  in  determining  whether  the  cloud  is  saturated  and  latent  heat  is  released. 
Therefore  this  will  be  an  area  for  further  research. 


4,0  Further  Work 

The  work  described  in  this  section  are  ideas  that  have  come  up  as  to  how  to  improve  the 
CRM  of  DELFIC.  Some  of  this  section  represents  ideas  proposed  by  other  authors.  Some 
of  this  .section  is  part  of  the  prospectus  presented  to  the  author's  research  committee  at 
AHT. 

4.1  Minor’s  recommendations 

It  was  stated  that  the  assumption  made  by  Norment,  that  the  specific  heat  of  entrained  air 
could  be  modeled  by  that  of  dry  air  (Norment,  1970),  may  not  be  valid.  This  would  require 
replacing  CPAI  with  a  new  variable  CPI  which  is  the  integral  of  Cp  instead  of  the  integral 
of  Cp,.  The  code  can  easily  be  altered,  but  the  way  Minor  defines  the  variable  CPI  is  in 
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FIGURE  4.  Effect  of  nmpHKt  wiad  jhear,  WIND  SHEAR  Comparisoo 


error,  The  new  vanable  should  be  CPI={CPAI+XE»CPWI)/(1.+XE)  not  CPI=(CPAI+X»- 
CPWI)/(1.+X)  as  stated  (Minor,  1988).  TTiis  new  variable  CPI  should  be  used  in  equations 
DERIV128  and  DERIV113  of  the  current  CRM. 

Miucr’s  other  area  of  concern  was  the  modeling  of  dry  and  wet  fallout  separately.  He 
states  that  P(t)  (or  F  in  Volume  I)  is  only  soil  and  not  total  fallout.  He  suggests  an  addi¬ 
tional  variable,  for  water  fallout.  He  stated  that  if  P(t)  included  both  dry  and  wet 

fallout,  the  particle  size  distribution  would  have  depended  on  water  content. 

In  reviewing  his  claims,  it  was  found  that  the  origmal  1967  DELFIC  did  account  for  con¬ 
densation  in  calculating  the  density  of  the  particles  and  their  volumes.  These  changes  in 
particle  de.nsity  and  volume  were  taken  out  with  the  first  revision  of  the  CRM  (Norment, 
1970).  In  addition,  the  ability  to  allow  condensed  wate.-  to  fall  out  of  the  cloud  was  elimi¬ 
nated  in  the  1979  version.  This  will  be  looked  at  in  conjunction  with  the  bigger  question  of 
fallout  in  general  mentioned  above. 

4.2  Heubsch’s  recommendations 

Heubsch's  1975  critique  of  the  revised  CRM,  stated  that  the  model  produced  physically 
unrealistic  results  or  results  that  disagreed  with  the  experimental  data.  Upon  closer  exami¬ 
nation  of  the  model,  Heubsch  proposed  the  errors  were  due  to  the  formulation  of  the  equa¬ 
tions  for  cloud  mass,  velocity,  temperature,  and  dimensions.  He  found  that  some  of  the 
terms  in  the  equations  contained  theoretical  errors. 

Heubsch  presented  an  amended  set  of  cloud  equatioits  which  removed  both  the  theoretical 
errors  and  output  diracpancies.  He  made  further  recommendations  on  how  to  improve  the 
CRM  and  proposed  three  levels  of  effort.  His  recommendations  are  important  since  he 
was  the  original  author  of  the  CRM.  However,  since  his  organization  no  longer  had  con¬ 
trol  of  DELPIC,  he  was  not  able  to  mandate  implementation  of  the  improvements. 

Norment  published  his  own  validation  study  of  the  revised  CRM  (Norment,  1977).  Nor¬ 
ment  did  change  the  momentum  (or  velocity)  equation,  as  Heubsch  had  recommended,  but 
left  the  other  equations  alone.  He  did  adjust  some  of  his  empirical  parameters,  simply  as 
another  method  of  reaching  better  agreement  with  experimental  data.  He  pointed  out 
which  parameters,  using  his  set  of  equations,  most  influenced  changes  in  the  results. 

Norment  did  not,  however,  change  the  other  equations,  in  particular  the  entrainment  equa¬ 
tion.  Heubsch  recommended  keqting  only  the  first  term  of  the  entrainment  equation  based 
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on  current  theory.  A  model  which  doe«  not  contrtdict  current  theory  and  matches  the  test 
date  is  needed. 


43  McGahan’s  recommendations 

In  addition  to  the  recommendation  to  go  back  and  coirect  the  violations  of  theory  in  the 
current  version  of  DELFIC's  CRM,  there  exists  an  outside  recommendation.  McGahan 
from  SAIC  is  currently  the  most  knowledgable  person  with  regards  to  the  theory  and  use 
of  the  DELFIC  model.  He  does  need  someone  to  review  the  CRM  for  its  errors,  but  in 
addidon,  he  would  like  to  improve  the  model.  He  wants  to  see  if  it  could  distribute  the  par- 
dcles  in  a  more  robust  method  than  the  ID  (verdcal)  placement  of  pardcles  it  now  uses. 
This  would  then  be  followed  by  a  comparison  with  the  higher  level  hydrodynamics  codes 
that  are  starting  to  be  used  for  nuclear  burst  modeling. 

His  recommendadon  involves  using  vortex  theory  and  a  two  dimensional  representadon 
of  the  cloud  pardcle  distribudon.  In  a  similar  fashion  to  how  DELFIC  distributes  the  pard¬ 
cles  in  ID  after  its  parcel  rise  subroutine,  vortex  theory  could  be  used  to  provide  a  2D 
flow-field  in  the  rising  cloud.  By  tracing  the  flow  of  pardcles  during  rise,  using  the  ou^ut 
values  from  the  cloud  rise  equadons,  the  posidon,  velocity,  and  acceleradon  of  particles 
could  be  tracked  up  to  stabilization. 

This  would  be  a  new  method  of  calculating  the  data  used  by  the  nuclear  effects  community 
which  is  cutrendy  available  only  from  the  supercomputer  hydrodynanucs  codes.  The  goal 
is  to  refine  the  verdcal  distribudon  of  the  larger  pardcles  in  DELFIC  results,  which  cur- 
rendy  stabilize  at  much  lower  aldtudes  than  the  hydrocodes  predict.  McGahan  thinks  that 
accounting  for  the  spadally  varying  flow  field  will  allow  a  better  treatment  of  the  large 
pardcle  tracking,  and  possibly  show  them  attaining  a  higher  aldtude.  The  implementadon 
of  this  would  also  provide  a  radial  distribudon  of  pardcles  that  is  not  cutrendy  available 
from  DELFIC.  Not  only  will  this  tracking  give  a  better  definidon  of  the  particle  distribu¬ 
tion  with  time,  but  also  give  dynamic  representadon  of  the  particles. 

4.4  DICE/MAZ  -  TASS  results 

/  hydrodynamic  codes  are  being  develc^>ed  more  fully,  their  output  is  being  trusted 
more.  Tne  idea  was  suggested  diat  the  output  of  TASS,  a  2-D  hydrocode,  may  want  to  be 
used  to  help  determine  the  best  values  for  the  remaining  parameters  in  CRM.  This  would 
be  in  addidon  to  the  validation  efforts  with  atmospheric  test  daut  This  area  is  still  being 
consideted  but  may  not  be  feasible. 
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5.0  Conclusion 


This  repo.t  showed  most  of  the  findings  during  a  recent  study  of  the  1979  version  of 
DELFIC’t  Cloud  Rise  Module  (CRM).  The  areas  cf  the  code  and  documentation  that  still 
contain  errors  were  listed  as  well  as  the  suggested  corrections.  Other  areas  of  possible 
improvement  to  the  CRM  found  using  a  higher  level  language  were  discussed.  Finally,  the 
proposed  work  for  the  current  research  at  AFTT  was  discussed.  These  areas  include  revis¬ 
ing  the  equations  in  the  CRM  to  be  self  consistent  with  conservation  laws  and  available 
data.  Also  the  inclusion  of  a  2-D  vortex  flow  field  for  particle  placement  within  the  cloud 
will  be  accomplished. 
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